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^ . ABSTRACT 

Galactic black hole binaries produce powerful outflows with emit over almost the 
entire electromagnetic spectrum. Here, we report the first detection with the Herschel 
observatory of a variable far-infrared source associated with the compact jets of the 
black hole transient GX 339—4 during the decay of its recent 2010-2011 outburst, after 
the transition to the hard state. We also outline the results of very sensitive radio 
ly-^ ■ observations conducted with the Australia Telescope Compact Array, along with a 

iy) ' series of near-infrared, optical (OIR) and X-ray observations, allowing for the first time 

£SJ , the re-ignition of the compact jets to be observed over a wide range of wavelengths. 

The compact jets first turn on at radio frequencies with an optically thin spectrum 
that later evolves to optically thick synchrotron emission. An OIR rcflare is observed 
. about ten days after the onset of radio and hard X-ray emission, likely reflecting 

the necessary time to build up enough density, as well as to have acceleration (e.g. 
through shocks) along an extended region in the jets. The Herschel measurements are 
consistent with an extrapolation of the radio inverted power-law spectrum, but they 
' highlight a more complex radio to OIR spectral energy distribution for the jets. 

, Key words: accretion, accretion disks - black holes physics - X-rays: binaries - ISM: 

jets and outflows - stars: individual (GX 339-4). 
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1 INTRODUCTION 

During outbursts of X-ray binaries, an enormous amount 
of ac cretion energ y is harnessed to produce powerful jets 
(e.g. iFended 120061 ). For a very long time, those jets were 
only detected i n rad io (e.g. lHiellming. Gibson fc Owenll975l ; 
iMirabel etHI Il992l ). However, the recent decade demon- 
strated that relativistic jets from black holes (hereafter BH) 
and neutron stars are mu lti-wavelength emitt e rs with emis- 
sion e ven up to 7 rays (jLaurent et alj l201ll : Corbel et al.l 
l2012bl ). Powerful self-absorbed compact jets (|FendeiH20oif ) 
are observed in the hard state of BH binaries, and are rem- 
iniscent of those found in active galatic nuclei. Models of 
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comp act jets (|Blandford fc Koniglll 19791 : iFalcke fc Biermannl 
Il995l ) predict that their spectral energy distributions (SED) 
should be a flat or slightly inverted power-law between 
the radio and a turnover frequency esti mated in the mid 
to near-infrared (NIR) fo r Galactic BHs dCorbel fc Fenderl 
120021 : iGandhi et all l201ll : iRussell et all l2012bh . Such NIR 
synchrotron excess emission from the jets has now been 
detected in a growing number of BH transients (e.g . 
Jain et al.ll200ll; iBuxton fc Bailvnll2004l: iBuxton et al.|[20ia: 



Dincer et aill2012l : IRussell et all 1201281 ). However, there is 



almost an unexplored four decades frequency gap between 
the radio and the NIR domains. Only very few millimeter 



and mid-infrared (IR) observations (e. g . I Fender et alj 



Migliari et ail 120071; iGallo et all 120071 : iRahoui et al 



2001 



2011 



Gandhi et alfellf ) have been conducted in this range and 
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usually with very limited sensitivity. New theoretical calcu- 
lations high light the possibility of a more complicated spec- 
tral shape (|Pe'er fc Case lla 2009), especially in the far to 
NIR range. Strong variability can be observed in the NIR 
on short timescale (~ seconds) allowing one to probe phys- 
ical changes at the base of the jet. These variations seem to 
vanish on longer timescales (above a few hundreds seconds , 
ICasella et aHl20ld : iGandhi et aHl201ll ; iRahoui et alll2012h . 

This points to the importance of sensitive sub- 
millimeter to far-IR observations that can only bridge the 
gap between the radio and the NIR ranges. In this Letter, 
we report the first far-IR observations with the ESA Her- 
schel satellite of the recurrent binary GX 339—4 during the 
decay of its recent outburst in 2011. GX 339—4 is believed 
to harbour a BH with a mass > 7 Mp, (jHvnes et al.ll2003l ; 



M unoz-Darias. Casares fe Martfnez-Paisl l2008h in a system 
with a low mass star companion, and located at a likely dis- 
tance of 8 kpc l|Zdziarski et al.ll2004l ). In Section 2, we de- 
scribe the Herschel observations along with complementary 
radio, NIR, optical and X-ray observations. The detailed 
multi-wavelength lightcurves, as well as the far-IR detec- 
tions, are also presented and these results are discussed in 
Section 3. 



2 OBSERVATIONS AND RESULTS 

2.1 Herschel far-infrared Observations 

Following the transition of GX 339—4 back to the hard state 
during the decay of its 2010-11 outburst, we triggered two 
DDT (Director Di scretionary Time) o bservations with the 
Herschel satellit^ (|Pilbratt et al.ll201Ch. GX 339-4 was ob- 
served with the PACS instrument ( Poglitsch et aljboich at 
70 /im and 160 fj,m on 2011, February 25 and March 6 (Ta- 
ble 0. We used a sequence of two so-called "mini-scan map" 
mode, providing an area of 50" in diameter with homoge- 
nous coverage for a total on-source time of 1152 s. The data 
were processed with HCSS version 6.0.1106 and we recon- 
structed the maps with the Tamasis software (Chanial et 
al. in prcpjfl As GX 339—4 sits in the Galactic plane in the 



Obs. ID 


Date 


f v (70 iim) 
(mjy) 


f v (160 fim) 

( mJy) 


1342215720 
1342215721 


2011-02-25 
2011-02-25 


21.0 ± 1.0 


13.2 ± 4.6 


1342216156 
1342216157 


2011-03-06 
2011-03-06 


13.9 ± 1.0 



1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 

2 see http://pchanial.github.com/tamasis-pacs/ 



Table 1. Herschel observations log of GX 339—4 

midst of a strong cirrus background, we paid a special atten- 
tion in the map reconstruction to preserve the structure of 
the background. The resulting flux density maps at 70 /im 
are shown in Fig. [T](a, b). A point source is clearly detected 
at the center of each map, as well as some extended emis- 
sion across the field. A single point source is visible in the 
difference image (Fig. [TJj) at the center, clearly indicating 
the detection at 70 fj,m of a variable source at the location 
of GX 339—4. All the other sources and structure have been 
subtracted out. In order to choose the optimum apertures 
for source and sky extraction, calibrate our photometry and 
check the impact of the varying cirrus background, we have 
run a set of simulations, where point sources were added to 
the bolometer timelines by applying our obser vation mode l 
to the observed PACS point spread function of lLutj (|2012f ) . 
We reconstructed the maps and performed aperture photom- 
etry on the simulated sources and GX 339—4, with various 
aperture size and sky annuli. We selected the combination 
giving the highest signal to noise for point sources (Table [1]). 
The cirrus confusion noise adds 0.77 mjy at 70 /im and 4.5 
mjy at 160 fim in quadrature to the point source sensitivity. 
As a result, we only obtain a 2.9 a detection at 160 fim by 
combining the two observations, while GX 339—4 is detected 
at high significance at 70 /im in each epoch. 

2.2 Radio, OIR and X-ray observations 

In order to probe the evolution of the SED from the compact 
jets, the Herschel observations were included in a large cam- 
paign of multi-wavelength observations. All radio observa- 
tions were performed with the Australia Telescope Compact 
Array (ATCA) located in Australia. The ATCA observa- 
tions were conducted at 5.5 and 9 GHz (also at 2.2 GHz dur- 
ing the first Herschel run) with the upgraded and sensitive 
CABB back-end (see Table 2 available as online support- 
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Figure 2. (Top): The lightcurves (with multiplicative factors as 
indicated in the legend) of GX 339—4 during the decay of the 
2010-11 outburst in radio (green), NIR (red), optical (purple), 
soft X-ray (filled black circle) and hard X-ray (open black circle). 
The arrows highlight the Herschel observations (observations # 
3 and 4). (Bottom): Evolution of the radio (green hexagon) 5.5 
to 9 GHz spectral index (see text) once GX 339—4 initiated its 
transition back to the hard state. The radio to far-IR spectral 
indices using the Herschel data) are represented by the red circle. 
The dashe d line indicates the transition back to the hard state ac- 
cording to lDincer et al. N2012I) . The vertical dotted line highlights 
the onset of OIR emission close to the peak of hard X-ray and 
radio emission. The numbers refer to the individual observations 
displayed in Fig. \3\ 



ing information) . The PCA and XRT detectors onboard the 
RXTE and Swift satellites provide the X-ray coverage dur- 
ing the soft to hard state transition. Further details on the 
radio a nd X-ray data analysis can be found in lCorbel et alj 
(1201231 ). The optical and near-infrared (herafter OIR) ob- 
servations were taken with the ANDICAM camera on the 
SM ARTS 1.3 m telesco pe located at Cerro Tololo in Chile 
fsee lBuxton et al"1l2012h . Images in H, J, I and V filters are 
obtained on a near-daily basis when GX 339—4 is visible in 
the night sky. The magnitudes have been converted to flux 
units using a total extinction of Av = 3.7 mag. 



2.3 The re-ignition of the compact jets 

Figure [2] shows the multi-wavelength lightcurves of 
GX 339-4 during the decay of its 2010-11 outburst. The 
dashed line indicates the increase of the hard X-ray emission 
on MJD 55594 associated with the transition of GX 339—4 
back to the hard state (see lDincer et al1l2012l ). Radio emis- 
sion from the com pact jets is quenched in the soft state 
|Fender et alJll999T ). Due to the increased sensitivity of the 
ATCA array and denser coverage, we fully observe, for the 
first time, the recovery of the radio emission during the soft 
to hard state transition when the compact jets are building 
up. Despite sparse initial coverage, the onset of radio emis- 
sion from GX 339—4 is clearly observed. It occurs in con- 
junction with the increase of hard X-ray emission (Fig. [2]). 
The peak of radio emission is concomitant with the peak in 
hard X-rays, as well as with the rise of the OIR emission. 
However, the overall shape of the reflares are different with, 



e.g., the hard X-rays appearing to be more symmetric than 
the radio peak. The peak flux in OIR is reached ~ 2 weeks 
after the maximum of radio emission. The timescale to reach 
the peak emission (assuming a transition on MJD 55594) is 
much longer (~ 17 days) in OIR compared to hard X-rays 
(~ 10 days), and maybe also compared to radio. We note 
that the transition may have even started in radio and X-ray 
as early as MJD 55578, if we consider the weak radio detec- 
tion (see Table 2 online) that is co nsistent with the le veling 
off of the X-ray timing properties jDincer et al.ll2012T ). 

If we define the radio flux density, F„, as F„ oc v a 
(with a the spectral index), we note a clear change, after 
the transition to the hard state (Fig. [2}, from a negative 
spectral index to a saturation around ~ +0.4. This trend 
corresponds to an evolution from optically thin to optically 
thick synchrotron emission, associated with the building up 
of the compact jets during the state transition. Connecting 
the radio measurements to the two Herschel observations at 
70 ^im allows us to precisely characterise the radio to far- 
IR spectral indices to 0.32 ± 0.01 on Feb. 25th and 0.38 
± 0.01 on Mar 6th, which is fully consistent with the radio 
estimates (Fig. [2]). The detections by Herschel of a variable 
source at the location of GX 339—4 during a hard state fur- 
ther indicate that this far-IR emission is likely related to the 
self-absorbed compact jets of GX 339—4. 



3 DISCUSSION 

Our campaign during the decay of the 2010-11 outburst of 
GX 339—4 allowed us to obtain for the first time a significant 
detection of the compact jets in the far-IR with Herschel, 
as well as to observe the full recovery of the compact jets 
in radio and OIR. This enables us to bridge the 4-decade 
frequency gap between the radio and OIR domains. 



3.1 Formation of the compact jets 

Inspection of the evolution of the broadband spectra (Fig. 
[3] ) demonstrates that the compact jets are first detected 
at radio frequencies as optically thin synchrotron emission 
(# 1 and 2). The radio spectrum becomes flat (Fig. [5] (bot- 
tom)) after the onset of the OIR reflare and its spectral 
index then increases up to a value of the order of +0.40, 
at a time when the OIR peaks and starts to decay. Such 
sequence has not been documented be fore in as much de - 
tails (e.g. the 1999 decay of GX 339-4. ICorbel et alj|200(ih . 
especially in light of the evolution of the radio emission 
during the soft to hard state transition. However, the pre- 
vious st udies that discuss some radio and/or OIR obser- 



vations ( Kalemci et al. 2005, 2006: Brocksopp et al.l 



Fender. Homan fc Belkmil 120091; IMiller- Jones et all 



2005; 



Iml; 



Dincer et al.ll2012l : IMiller- Jones et al.ll2012T ) are consistent 
with our new findings. The comparative evolution of the 
radio and OIR lightcurves clearly points to a delay (~ 2 
week s) between these two domains (see also ICoriat et al.l 
l2009f ). In fact, this is consistent with the reverse trend that 
is observed during the hard to soft state transition, with an 
initial quenching (by one to two weeks) of the NIR emission 
jHoman et al l 120051: iRussell et al.ll2008l : ICoriat et all 120091 : 
iBuxton et al.ll2012l : lYan fc Yull2012r i~ 

In the model of self-absorbed compact jets 
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(|Blandford fc Konigll 1 19791 ; iFalcke fc Biermannl ll995T) 



adap ted for stellar mass BHs (e.g. lMarkoff, Nowak fc Wilms! 
l2005[ l. the OIR emission arises from closer to the base of 
the jets, whereas the radio emission comes further down the 
jets. The observed delay between OIR and radio emission 
during the formation is not consistent with any travel 
time within the jets. We still do not fully understand the 
internal processes in jets leading to particle acceleration 
and ultimately radiation. However, this delay is likely 
associated with the timescale for building up turbulences 
and strong shocks along an extended zone in the magnetized 
jet plasma, as well as for reaching high enough densities for 
the flow to become optically thick at OIR frequencies. One 
can imagine the following scenario to explain the evolution 
of the broadband spectra. Once the jets are launched, an 
initial zone of particle acceleration, associated with the 
first shocks, forms. While the density of the jet plasma 
is still low, this results in an optically thin synchrotron 
power-law spectrum peaking in the radio band. As the 
density in the jets increases, this leads to a transition to 
higher optical depth. At the same time, pile-up of the 
flow or internal shocks (|Malzad l20l3) build up and cause 
particle acceleration to occur closer to the base of the jets. 
Alternatively, an increase in particle acceleration efficiency 
at the jet base could also boost the energy of the particles, 
giving rise gradually to enhanced synchrotron emission up 
to OIR. 

One would therefore expect the observed behaviour: a 
transition from a steep to a flat spectrum with an evolution 
of the peak power from radio to NIR. A reverse scenario (de- 
crease in the particle densities and/or the efficiency of parti- 
cle acceleration) would also be consistent with the trend (ini- 
tial OIR quenching) observed during the hard to soft state 
transition. The onset of radio emission can thus be seen as 
a tracer (better than OIR) of the jets re-ignition, while the 
broadband evolution leading to inverted synchrotron spec- 
trum peaking in the OIR traces the development of the jets 
power and their internal structure during outbursts. 

3.2 On the far- to near-infrared emission 

The very good sensitivity of the Herschel/PACS instru- 
ment and of the upgraded ATCA demonstrates that the 
far-IR emission of GX 339—4 is clearly and significantly 
above a simple extrapolation of the radio to NIR spec- 
trum (see broadband spectral evolution in Fig. during 
the phase of compact jets formation). This is in contradic- 
tion with the usual assumption of a simple powerlaw (e.g. 
iMarkoff. Nowak fc W ilms 2005]) with a turnover frequency 
in the mid to near-IR range. Subtracting the thermal compo- 
nent in the OIR would further increase the far-IR emission 
relative to the OIR. This is in fact reminiscent of the 850 
Ami (350 GHz) JCM T detection of XTE J1118+480 in 2000 
ijFender et al.ll200ll ). which was also clearly above the radio 
to OIR inverted spectrum (see their Fig. 2). The radio to 
sub-millimeter spectral index of +0.5 for XTE J1118+480 
is furthermore quite comparable with the value measured 
during our Herschel observations of GX 339—4 (Fig. [2]). 
V404 C yg may display occ asionally a similar broadband 
spectra (|Russell et alj|2012bl ). 

These three sources seem to point out a more com- 
plex broadband spectra for the compact jets during some 
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Figure 3. Evolution of the radio to OIR broadband spectra dur- 
ing the soft to hard state transition. The numbers (and associated 
colors) refer to the epochs highlighted in Fig. [2] The dashed lines 
corresponds to the extrapolation of the radio spectra (with dotted 
lines to indicate the 1-<t error when it is not too large). For the 
spectra # 3 and 4 (which includes the Herschel observations), 
we use the measured radio to far-IR spectral index. The open 
blue circle is the 160 /an detection combining the two Herschel 
observations. Inset: Zoom of the OIR spectra # 1, 3 and 6. The 
dashed lines correspond to an extrapolation of the V-band w ith 
an index of +1.5 (possible for an irradiated disk, lHvnes!l2005l l for 
a selection of the spectra. 



epochs. The optically thick synchrotron emission from the 
jets manifests itself by an inverted power-law from the ra- 
dio up to at least the far-IR (same spectral index, Fig. [2}. 
The turnover frequency during our two Herschel observa- 
tions is likely above 5 x 10 12 Hz, which is consistent with 
previous estimates. It is cle ar from the evoluti on of the OIR 
lightcurves (Fig. [2] and also IDincer et al1l2012t l that most of 
the NIR emission during the refiare is not from the intrinsic 
decaying thermal emission from the accretion disk. Further- 
more, this OIR refiare can not be completely due to the 
optically thin synchrotron emission from the compact jets, 
because it displays a rather flat spect ra (after subtractio n 
of the thermal emission) according to IDincer et al.l l)2012h . 
This means that an other emission process may dominate th e 
OIR range (see also lRahoui et alj|2012l ; IDincer et alj|2012h . 
Even if the OIR refiare is delayed compared to the radio syn- 
chrotron emission, this OIR emission could still be related 
to the jets via, e.g., pre-shock syn chrotron from the thermal 
parti cles at the base of the jets ()Markoff. Nowak fc Wilmsl 
l2005h . The OIR could also be affected by an episode of parti- 
cle ac celeration in the jets i n a particularly strong magnetic 
field l)Pe'er fc Casellal [20091 1. Alternatively, the OIR excess 
may be re lated to synchrotron emission from the hot accre- 
tion flow (|Veledina. Poutanen fc Vurm|l20l"lh . At first sight, 
the shape of the OIR spectra (see inset in Fig. [3} could 
also possibly highlight, as a viable alternative mechanism, 
the role of irradiation of the disk (Hynes 2005) by the jets, 
which form during this period. 

With optically thick synchrotron emission in the 10 9 - 
5 x 10 12 Hz (~ 70 /im) range and assuming a radiat i ve eff i- 
ciency of e= 0.05 [see discuss i on inlCorbel fc Fenderl l)2002h . 
or recently in lNemmen et all (|2012h who suggest that e could 
be as high as 0.15], we obtain a likely lower limit on the 



power of compact jets during Herschel observations #1 of 
~ 10 35 erg s _1 . Adding an optically thin jet component (of 
spectral index -0.6) in the 5 x 10 12 to 10 14 Hz increases the 
jet power to ~ 10 36 erg s _1 . This corresponds to at least 
17 % of the 0.1-200 keV luminosity (6 x 10 36 erg s" 1 for 
a distance of 8 kpc), a significant fraction of the accretion 
power (assuming no correction for bulk relativistic motion). 
This probably implies that the jet spectral break can not 
be much higher than the frequency of the Herschel observa- 
tions. For a comparison, in the brighter hard state during 
the rise (with an X -ray luminosity of 1.5 xlO 38 erg s _1 ), 
I Gandhi et all l|201ll ) measured this break at ~ 5 x 10 
Hz. A jet frequency break moving to lower frequency is ex- 
pected as the mass accretion is r educed during the outburst 
decay (iHeinz fe Sunvaevl 120031 : iFalcke. Kording fe Markofj 
120041 ; iRussell et al.ll2oT2bT ). Furthermore, the Herschel mea- 
surements would translate into a magnetic field strength of 
~ 2 x 10 3 G in the acce leration zone and a jet cross sec - 
tion radius of ~ 10 10 cm (|Chatv. Dubus fe Raichoorll201l[ ). 
These numbers would eventually imply compact jets with 
weaker magnetic field and a larger base during the decaying 
hard state compare d to the rising hard state observed by 
I Gandhi et all (|201ll ) (if we take the same assumptions). 

To summarize, our new campaign highlights a much 
more complicated radio to OIR spectral shape than usually 
assumed, especially around the time of the transition. How- 
ever, on many occasions, the OIR measurement s may still lie 
on the extrapolation of th e radio spectrum (e.g. lCoriat et all 
120091 ; iBuxton et alj|20l3 ). This reinforces the need for new 
observations, e.g. with the ALMA array or with sensitive 
mid to far IR observations, to bridge the gap between the 
radio and the OIR domains and constrain the evolution of 
the jets SED and power during BH outbursts. 
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